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■ We show that it is now possible to fully determine the CKM matrix, for the first time, 

O . using lattice QCD. \Vcd\-, iV'csl, \Vub\-, \ycb\ and \Vus\ are, respectively, directly determined 

Q i with our lattice results for form factors of semileptonic D — > ttIv, D Klv, B irlv, 



B — > Dlu and K ttIu decays. The error from the quenched approximation is removed 
by using the MILC unquenched lattice gauge configurations, where the effect of n, d and 
s quarks is included. The error from the "chiral" extrapolation (mi rriud) is greatly 
reduced by using improved staggered quarks. The accuracy is comparable to that of the 
Particle Data Group averages. In addition, \Vud\, \Vtb\, \Vts\ and \Vtd\ are determined by 
, using unitarity of the CKM matrix and the experimental result for sin (2/3) . In this way, 

I we obtain all 9 CKM matrix elements, where the only theoretical input is lattice QCD. 

We also obtain all the Wolfenstein parameters, for the first time, using lattice QCD. 
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rin : 1 Introduction 

^ • The Cabibbo-Kobayashi-Maskawa (CKM) matrix, which relates the mass eigenstates and the 
weak eigenstates in the Standard Model electroweak theory, is a set of parameters. To deter- 
mine each CKM matrix element, one requires both theoretical and experimental inputs. On the 
^ ' theoretical side, one needs to know relevant hadronic amplitudes, which often contain nonper- 
c3 ■ turbative QCD effects. A major role of lattice QCD is to calculate such hadronic amplitudes 
reliably and accurately, from first principles. One can then extract the CKM matrix elements 
by combining lattice QCD as the theoretical input with the experimental input such as decay 
rates. In this paper, we show that it is now possible to fully determine the CKM matrix, for 
the first time, using lattice QCD. The result for the full CKM matrix with lattice QCD is: 

/ \Vud\ \Vus\ \Vub\ \ 

0.9744(5)(3) 0.225(2)(1) 3.5(5)(5) x IQ-^ 

\Vcd\ \Vcs\ \Vcb\ 

0.24(3)(2) 0.97(10)(2) 3.9(1)(3) x lO^^ 

\Vtd\ \Vts\ \Vtb\ 

V 8.1(2.7) X 10-3 3.8(4)(3)xl0-2 0.9992(0)(1) / 

where the first errors are from lattice calculations and the second are experimental, except the 
one for \ Vtd\ which is a combined lattice and experimental error. The results for the Wolfenstein 
parameters with lattice QCD are: 

A= 0.225(2)(1), A = 0.77(2)(7), p = 0.16(28), r/ = 0.36(11). (2) 
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To directly determine 5 CKM matrix elements |Ks|, \Vub\, \Vcb\ and |V^s|), we use 

5 semileptonic decays {D ttIu, D Klu, B ttIu, B Dlv and K ttIu), for which 
the techniques for lattice calculations are well established, and thus reliable calculations are 
possible. The accuracy of previous lattice calculations was limited by two large systematic 
uncertainties — the error from the "quenched" approximation (neglect of virtual quark loop 
effects) and the error from the "chiral" extrapolation in light quark mass {mi mud) ■ Both led 
to effects of around 10-20%. Our present work [H |2j successfully reduces these two dominant 
uncertainties. The error from the quenched approximation is removed by using the MILC un- 
quenched gauge configurations [3 , where the effect of u, d and s quarks is included (nj = 2 + 1). 
The error from the chiral extrapolation is greatly reduced by using improved staggered quarks. 
With this improved approach, the accuracy of the 5 CKM matrix elements is comparable to 
that of the Particle Data Group jlj. The results for \Vub\, \Vcb\ and \Vus\ are preliminary. 

Since we determine all 3 elements in the charm row of the CKM matrix, \Vcq\ {q = d, s, b), 
we can check a unitarity condition on the CKM matrix using only our results from lattice QCD. 
Adding in \Vus\ and \Vub\, we then use CKM unitarity to determine the other 4 CKM matrix 
elements (iKidl; \Vtb\, \Vts\ and |l^d|)- In this way, we obtain all 9 CKM matrix elements and 
all the Wolfenstein parameters, where the only theoretical input is lattice QCD.^ 

This paper is organized as follows. In Sec. El we present our results for the 5 semileptonic 
decays and the 5 CKM matrix elements. In Sec. IHlwe first check unitarity of the charm row of 
the CKM matrix, and then give the results for the other 4 CKM matrix elements from unitarity, 
as well as the Wolfenstein parameters. In Sec. ID we give a summary and discuss future plans. 
The results for the D and B decays have been presented in Refs. IHHI. This work is a part of 
our ongoing project of heavy quark physics with lattice QCD; see also j3ini- 



2 5 CKM matrix elements from 5 semileptonic decays 

2.1 D tt{K)Iv, \V,d{s)\ and B ttIv, \Vub\ 

The differential decay rate dV /dq^ for the heavy-to-light semileptonic decay H —>■ Plv is pro- 
portional to |V^jp|/-(.(g^)p, where /+ is a form factor of the relevant hadronic amplitude defined 
through 

{P\V^\H) = U{q'){pH+pp-Ar + fo{q')A^ (3) 



= [v^f\\{E)+plf^{E)_ . 

Here q = Ph — Pp, A.^ = (jnu ~ "^p) q^ 1^ i = PhI'^h-, P± = Pp — Ev and E = Ep. To 

2 

determine the CKM matrix element \Vij\ with the experimental rate J'^r'"'' dq^ {dT/dq'^), we 
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calculate /+_o as a function of q^. Below we briefly describe our analysis procedure 

We first extract the form factors f\\ and f±, and carry out the chiral extrapolation in mi for 
them at fixed final-state energy E. To this end, we interpolate and extrapolate the results for 
/|l and f± to common values of E using the parametrization of Becirevic and Kaidalov (BK) 
[7j. We perform the chiral extrapolation [mi —>■ mud) at each E using the NLO correction in 

"'^We assume here that the Standard Model is correct. We plan to test the Standard Model in future work. 
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Figure 1: Form factors from BK-based (filled symbols and curves) and non-BK-based (open 
symbols) analyses for D ttIu (left panel) and B — >■ ttIu (right) decays. 

staggered cliiral perturbation theory |B^. We try various fit forms for the extrapolation, and 
the differences between the fits, 2~4%, are taken as associated systematic errors. 

We then convert the results for f± and f\\ at m/ = rriud, to /+ and /q. To extend /+ and /o 
to functions of q^, we again make a fit using the BK parameterization 0, 



(1-^(1 
where (f = q'^/mjj,. We obtain pj|(2] 

= 0.64(3), a 



for the D decay, and 



+ 



0.73(3), 



0.23(2), 



a 



0.44(4), 
: 0.50(4), 
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0.63(5), (3 



B^TT 



: 1.41(6), 

= 1.31(7), 
1.18(5), 



(4) 

(5) 
(6) 

(7) 



for the B decays, where the errors are statistical only. To estimate the error from the BK 
parameterization, i.e., the error for dependence, we also make an alternative analysis, where 
we perform a 2-dimensional polynomial fit in {mi, E{q^)) . A comparison between two analyses 
is shown in Fig. ^ The results for D decays agree well with recent experimental results p. 

We now determine the CKM matrix elements by integrating \ f+{q^)\^ over q^ and using 
experimental decay rates jHHnHH]. For \Vub\ we use a combined average of the decay rate for 
g2 > 16 GeV2 in Refs. Hm and [H]. We obtain 



\VJ\ = 0.239(10)(24)(20) , l^sl = 0.969(39)(94)(24) 



from the D decay, and 



IV;^! X 10^ = 3.48(29) (38) (47) 



(9) 



from the B decays, where the first errors are statistical, the second systematic, and the third 
are experimental errors from the decay rates. The systematic errors are dominated by the 
finite lattice spacing effects, i.e., the lattice discretization effects; see Table [T] The results for 
the CKM matrix elements agree with the Particle Data Group averages |3] with a comparable 
accuracy. A similar calculation by the HPQCD collaboration is also underway 
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2.2 B Dlu, \Vcb\ and K ttIu, |Ks| 

The form factors oi B —>■ Dlu and K —>■ -n-lu decays can be calculated more accurately than 
those of heavy-to-light {H Plv) decays, due to the symmetry between the initial and final 
states. The differential decay rate of i? ^ Dlv is proportional to the square of \Vch\J^{w), where 
T{w) is a linear combination of hj^{w) and h-{w) defined through 

{D\V^\B) = y/mBruD x [h+{w){v + vY + h_{w){v - v')"], (10) 

with V = ps/mB, v' = Pd/^d and w = v ■ v' . To extract iKfel, we calculate the form factor at 
w = 1, ^(1), by employing the double ratio method ^31- The light quark mass dependence of 
^(1) is mild, and by extrapolating the result linearly to — we obtain [2] 

J^B^'^^l) = 1.074(18)(16), (11) 

where the first error is statistical, and the second is systematic, as summarized in Table ^ 
Using our result Eq. (jllj) and an average of experimental results for |Vc{,|J^(l) ^14j, we obtain 

iKbl X 10^ = 3.91(07)(06)(34), (12) 

where the errors from the lattice calculation (first two) are much smaller than the experimental 
one (third). 

Finally we study the K tcIu decay to determine \Vus\- The expression for the K ^ n 
decay amplitude is given in an analogous way to Eq. Q. We calculate the K'^ — tt~ form 
factor at = 0, /+(0) = /o(0), by employing the following three steps method, as in Ref. P3]: 

1. Extract the /o form factor at = q'^^^ by applying the double ratio method similar to 
the 5 ^ D case 

2. Extrapolate /o at = q^^^ to = by using the experimental result [12] for the slope 
parameter Aq defined by /o(g^) = /o(0)/(l — Aog^).^ This gives /o(0) = /+(0) for each mi. 

^An unquenched lattice calculation of Aq should be done in the future. 



semileptonic decay 


D n{K)lu B - 


tt/z/ B - 


Dlu K - 


— > vr/z/ 


CKM matrix element 


Vcd{s) 


iKbl 


\Vcb\ 


\Vus\ 


g^ dependence 


2% 


4% 




<1% 


nil ^ud extrapolation 


3%(2%) 


4% 


1% 


1% 


operator matching 


<1% 


1% 


1% 


<1% 


discretization effects 


9% 


9% 


<1% 




total systematic error 


10% 


11% 


2% 


1% 


error in PDG average 


5%(1%) 


13% 


4% 


1% 



Table 1: Systematic errors in lattice calculations. For comparison, the error for each CKM 
matrix element by the Particle Data Group j3] is shown in the last row. 
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Figure 2: mi-dependence of A/ for K ttIu decay, together with results from Refs. [TH| I17j. 

3. Perform the m/ — rriud extrapolation for /+(0). To this end, we subtract the leading 
logarithmic correction /2 in chiral perturbation theory, i.e., define A/ = 1 + /2 — /+(0). 
We make a fit to A/ adopting an ansatz. A/ = {A + Bmi){ms — miY, where A, B are fit 
parameters. The (m^ —miY dependence is expected due to the AdemoUo-Gatto theorem. 

As an exploratory study of f^^^, we use an improved staggered action for u, s quarks and an 
improved Wilson action for the d quark. The mrdependence of A/ and the extrapolated result 
are shown in Fig. |21 together with a recent quenched lattice result ^3] and an earlier result by 
Leutwyler and Roos jUj. Our preliminary result is A/ = 0.015(6) (9), giving 

/f---(0) = 0.962(6)(9), (13) 

which agrees well with those of Refs. ITT] . Combining with a recent experimental result for 
|Kts|/+(0) [in], we obtain 

IK.I = 0.2250(14)(20)(12). (14) 



• Nf=2+1 

□ Nf=0 (Becirevic et al '04) 
o Leutwyler&Roos'84 




3 Other 4 CKM matrix elements using unitarity and 
Wolfenstein parameters 

Having the 5 CKM matrix elements directly determined from the 5 semileptonic decays, we 
can check unitarity of the second row of the CKM matrix. Using Eqs. (jH)) and (fT^ . we get 

(iKdP + iKsP + IK.P)'/' = 1.00(10)(2), (15) 

which is consistent with unitarity. Hereafter the first error is from the lattice calculation and 
the second is experimental, unless otherwise stated. 

We now use unitarity of the CKM matrix to determine the other 4 CKM matrix elements. 
\Vud\, \Vtb\ and \Vts\ are easily determined: 

\Vud\ = (l-lKsl'-lK;,!')'/' = 0.9744(5)(3), (16) 
\Vtb\ = (l-|K6|'-|K.fer)'/' = 0.9992(0)(1), (17) 
iVtsl = \V:,Vub + V:M / \Vth\ ~ K:V,t\/\Vtk\ = 3.8(4)(3)xl0-^ (18) 
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Before proceeding to extract \Vtd\, let us give some of the Wolfenstein parameters. Using 



Eqs. dUl), (O and ©, we get 

A = iKsl = 0.225(2)(1), (19) 

A = \V,,\/X' = 0.77(2)(7), (20) 

(p2 + ^2)1/2 ^ \V^,\/{AX') = 0.40(6)(6). (21) 



It is difficult to extract \Vtd\ accurately using unitarity alone, because Vtd (= AX^{l—p—iri) in 
Wolfenstein parameterization) contains a large CP-violating phase and its magnitude is small. 
We use here the experimental result for sin(2/3) from B — »• {cc)K^*^ decays. By performing a 
unitary triangle analysis with sin(2/3) = 0.726(37) jHj and Eq. (j2H), we obtain 

p = 0.16(28), r] = 0.36(11), (22) 

\Vtd\ = 8.1(2.7) X 10-3 (23) 

with a combined lattice and experimental error, completing the full CKM matrix. 

4 Summary and future plans 

In this paper, we have shown that it is indeed possible to fully determine the CKM matrix 
using lattice QCD. \Vcd\, \Vcs\, \Vub\, \Vcb\ and \Vus\ are directly determined with our lattice 
results for form factors of 5 semileptonic decays. The error from the quenched approximation is 
removed by using the MILC unquenched gauge configurations. The error from the m/ — >■ rriud 
extrapolation is significantly reduced by using improved staggered quarks. The accuracy is 
comparable to that of the Particle Data Group averages, as seen in Table H The other 4 CKM 
matrix elements \Vud\, \Vtb\, \Vts\ and \Vtd\ are then determined by using unitarity of the CKM 
matrix and the experimental result for sin (2/5). In this way, we obtain the full CKM matrix, 
Eq. (P), and the Wolfenstein parameters, Eq. We emphasize that the only theoretical input 
in this paper is lattice QCD. This concretely shows that lattice QCD is one of the most powerful 
tools in aiding our understanding of fiavor physics. 

As for future works, we are planning to improve the accuracy of form factors of D ^ 'k{K)Iv 
and B tcIu decays, which is currently 0(10%). To reduce the largest systematic error from 
the lattice discretization effects, we are repeating the calculations at a smaller lattice spacing. 
The unquenched calculations of leptonic decays (-D(s) — > li^, B(^s) — >• If) are also underway [HIITH]. 

To test the Standard Model, we need to calculate mixing parameters of neutral B and K 
mesons, Bb and Bk- The unquenched {nj = 2 + 1) calculation of Bb and Bk is underway by 
the HPQCD collaboration JH], and is planned by the Fermilab Lattice collaboration and by 
the LANL lattice group. The initial accuracy will be 0(10%), and we hope to have an accuracy 
of 5% in the near future. We will then be able to have tighter constraints on the (p, rj) plane, 
leading to a precision test of the Standard Model, and guiding the search for new physics. 
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